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Abstract. The first part of the paper discusses nuclear properties relevant to predict 
compound reactions. The second part addresses direct reactions with special emphasis 
on direct neutron capture. 



1. Introduction 

Reaction rates are central in the modeling of nucleosynthesis and astrophysical energy 
generation and are computed by an integration of reaction cross sections across an 
energy range defined by the plasma temperature. Due to the comparatively low energies 
(E < 10 MeV for charged projectiles and E < 1 MeV for neutrons) relevant in 
astrophysical plasmas, mainly three reaction mechanisms are important pQ : compound, 
resonant, and direct reactions. In principle, both resonant and direct mechanism may 
contribute to the reaction cross section although one of them is negligible in most 
cases. The nuclear level density in the compound nucleus, at the energy at which 
the compound nucleus would be formed, determines which mechanism is dominating [2]. 
The statistical Hauser-Feshbach model can be applied when there are many, overlapping 
resonances which can be treated with averaged transmission coefficients. In the absence 
of resonances and strong tails of resonances, the direct mechanism will be important 
even at low energy pQ. In between the two regimes, the contributions of single resonances 
require the precise knowledge of their properties, including energy, spin, parity, as well 
as the phases of wave functions to account for possible interference effects. 

Describing the interplay of different mechanisms has been a longstanding problem in 
the study of nuclear reactions and requires precise predictions of excited state properties. 
A thorough discussion is beyond the scope of this paper but it is emphasized that 
this, among other things, shows the great importance of measuring cross sections in 
the relevant energy range. Here, we are focusing on some of the ingredients necessary 
for calculations in the Hauser-Feshbach model. One of the main difficulties in the 
determination of reaction rates for astrophysics lies in the fact that most reactions 
involve unstable nuclei which are not (yet) accessible in the laboratory. Therefore, the 
required properties cannot be extracted from experiment but rather have to be predicted 
by theory. In the second part of the paper we discuss direct capture and compare the 
results to Hauser-Feshbach cross sections for selected nuclides. 
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2. Compound reactions 

2.1. Important nuclear properties 

Averaged widths (or transmission functions) are the central quantities in the calculation 
of cross sections in the Hauser-Feshbach model [3]. The relevance of the nuclear level 
density not only for the identification of the reaction mechanism but also for the 
computation of transmission functions has been discussed elsewhere [U EJ Hj. The 
sensitivity of the results to variations in the level density is not strong because i) 
transitions to low lying states are dominating and these are usually explicitly included 
[5. 6J, and ii) the thermal population of target states in the plasma allows for additional 
transitions washing out the effect of non-equally distributed parities [I]. 

Averaged radiative widths are usually calculated including El and Ml transitions 
[3 |5j ||, more recently also E2 j9]. For El transitions, Lorentzian shapes with 
modified low-energy tails [5] or strength functions from QRPA calculations [10J are 
used. Frequently, the appearance of low-lying additional GDR strength is discussed, 
e.g. stemming from soft-mode vibrations, often termed "pygmy resonance". This 
could, in principle, strongly enhance neutron capture cross sections of very neutron- 
rich nuclei [U [10]. However, it will only be important for r-process nucleosynthesis 
if three conditions are met: i) the additional strength has to be located below the 
neutron separation energy of the nucleus because otherwise there will be no 7 transitions 
affected, ii) the statistical model has to be applicable because otherwise direct capture 
may be dominating which does not see the GDR+pygmy strength, and Hi) the involved 
nuclei should not participate in reaction equilibria ((n,7)-(7,n) equilibrium) but rather 
be produced by individual captures or photodisintegrations. The last item calls for 
nuclei located outside the main r-process path which are reached in the brief freeze- 
out phase. There is no consistent picture yet. Previous microscopic calculations 
found some strength below the neutron threshold [TT], while recent calculations find 
the pygmy resonance consistently several MeV above the neutron separation energy 
[T2l [T3t HU [15j [16] . The latter would rule out an effect on astrophysical rates provided 
the width is not so large as to reach below the separation energy. The location of the 
additional strength above the separation energy for Sn isotopes is confirmed by a recent 
experiment [TTj, [181 EH]- Earlier experiments found additional El strength below the 
separation energy in semi- magic isotopes of Ba, Ce, and Sm [20]. Future experiments 
have to test the predictions further away from stability. 

Another complication in the prediction of 7 transitions arises due to isospin selection 
rules suppressing El and Ml T = — > T = transitions. This leads to a strong 
suppression of (0,7) cross sections of self-conjugate nuclei [21] and some suppression of 
(p,7) and (n,7) on such nuclei. Often, this effect is included by introducing arbitrary 
suppression factors [8j [22] because the original Hauser-Feshbach equation implicitly 
assumes complete isospin mixing [3] . It is possible to generalize the model to explicitly 
treat isospin mixing [3, [21] . 

Nuclear masses determine the reaction Q-values in the open channels. Beyond the 
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region of measured masses, theoretical models have to be employed. Special attention 
has to be paid to the transition from the measured to the unmeasured region when 
computing Q-values. Special care has to be taken to avoid artificial breaks in the 
relative channel energies when subtracting values from different sources. 

2.2. Optical potentials 

Transmission functions for particles are usually predicted by utilizing the optical model 
PQ. This requires an optical potential with real and imaginary parts that are, in 
principle, dependent on the type of the projectile, target mass, and projectile energy. 
Because of the sensitivity of the transmission coefficients to the depth and geometry of a 
potential, the calculated cross sections are also very sensitive to the choice of potential. 
Global potentials are necessary for large-scale predictions which should also give reliable 
results off stability where no experimental information is available. 

The largest uncertainties in this respect are in optical potentials for charged 
projectiles due to the low incident energies occurring in astrophysical applications. The 
relevant energies are close to the Coulomb barrier. This makes it extremely difficult 
not just to measure cross sections but also to determine optical potentials by scattering 
[23] . Available global potentials are usually derived at higher energies and/or only for 
a limited mass range. 

The microscopic optical potential of [21] (JLM) with low-energy modifications 
by [25] is widely used for neutrons and protons. With a few exceptions it has been 
successful in reproducing experimental reaction data across a large range of target 
masses. Recently, new parameterizations of the JLM potential have been derived [261127] 
with the aim to improve its performance at energies above 160 MeV. However, this new 
potential does not work at low energies, as has been found in several comparisons to 
recent data [2E1 [29], whereas the JLM potential performed much better. However, it 
has been argued that the isovector components of the JLM potential may be too weak 
[26] . Indeed, in a comparison to recent 70 Ge(p,7) and 76 Ge(p,n) data [29] we find that 
the absolute value and the energy-dependence of the S"-factors of this reaction and of 
previously measured (p,7) and (p,n) reactions on Se and Sr isotopes can be better 
reproduced when increasing the JLM imaginary depth by 70%. This can be attributed 
to a stronger isovector component in the imaginary part. For further details, see |29j . 

Regarding optical potentials for a particles, see [30] for an extended discussion. 
There are no microscopic potentials available and attempts to derive a global potential 
are restricted to fits of Woods-Saxon shapes. Interestingly, despite of the higher 
Coulomb barrier, a captures at low energies can still be described within a factor of two 
by the simple potential of [31] which was derived from a scattering at energies of several 
tens of MeV. However, a trend is found in measured values which often are well predicted 
theoretically closely above the Coulomb barrier height and exhibiting a slowly widening 
difference at sub-barrier energies (see, e.g., [321 133]). This indicates that the barrier 
transmission may not be described correctly and that an additional dependence on 
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barrier height may have to be introduced in the potentials. Another problem may arise 
due to the formation of a clusters in the compound nucleus. The difficulty in describing 
low-energy a-induced reactions is similar to the well-known problem of predicting half- 
lives of spontaneous a-emitters. 

There is one remarkable case, 144 Sm(a,7), showing the by far strongest discrepancy 
between experiment and theory at the lowest measured energy of 10.2 MeV [31]. Because 
of the importance of the extrapolation to 9 MeV for 7-process nucleosynthesis, it is 
urgently required that this reaction be remeasured at sub-Coulomb energies! 

In general, when performing measurements aiming at an improvement in optical 
potentials one has to select reactions properly. The potential of interest should give rise 
to the smallest width of the involved channels because otherwise the sensitivity of the 
cross section to a variation of the potential is small [29l [35] . For a potentials this can be 
achieved by studying (a, 7) around the Coulomb barrier or in (a,n) and (n,a) reactions. 
For the proton potential, (p,n) or (n,p) reactions can be useful. This has also to be 
considered when fitting optical potentials to reaction data. Only reactions sensitive to 
the optical potential will provide useful constraints. Such fits including reaction data 
may be of particular importance because it has been pointed out recently that there 
may be a fundamental difference between potentials derived from scattering on a "cool" 
target and emission from a "hot" nucleus [36]. A dependence on nuclear temperature 
may have to be included in the potentials. 

3. Direct capture 

3.1. Importance for neutron-rich nuclei 

The vast majority of reactions relevant for astrophysics can be treated in the statistical 
model. However, the nuclear level density is low in nuclei which are light (^4 < 20) or 
have closed shells. Towards the driplines, the neutron- or proton-separation energies 
are decreasing which leads to formation of the compound nucleus at low excitation 
energy and thus with a low level density at that formation energy. In all these cases, 
the Hauser-Feshbach model will not be applicable anymore and has to be supplemented 
by accounting for resonant and direct processes. This is straightforward provided the 
required nuclear properties are known and the mechanisms are clearly separated. The 
latter will be the case when one of the reaction mechanisms dominates. If two of the 
mechanisms exhibit similar strength, optical potentials may have to be modified to 
account for the increased flux into the other reaction channel. Additionally, coherent 
summation of the channels may be required to allow for interference effects. 

Even in heavy, closed shell nuclei direct capture may have a non-negligible 
contribution as, e.g., was demonstrated in [37] where a discrepancy between an activation 
and a resonance measurement of 208 Pb(n,7) was shown to be caused by this mechanism. 
It becomes dominant in the absence of resonances. More specifically, it will be non- 
negligible for targets close to the driplines. Such nuclides are produced in the rp-process 
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Figure 1. Relation between direct neutron capture and compound capture as function 
of neutron separation energy for isotopes of Ti, Cr, Fe [38], Sn [39], and Pb [37] 



on the proton-rich side and the r-process at the neutron-rich side of the nuclear chart. 
Figure Q] shows the increasing importance of direct neutron capture with decreasing 
neutron separation energies. (It is to be noted that these estimates just compare 
compound and direct cross section without accounting for interference effects.) 

3.2. Statistical direct capture 

It has been shown [37J that the direct cross section is very sensitive to the predicted 
properties of the final states and that different microscopic models yield vastly different 
results. The sensitivity is much higher than in the case of compound reactions because no 
averaged quantities are used. In order to circumvent the problem of the exact prediction 
of states it has been suggested [HH HTJ H2] to employ averaged properties also for direct 
captures, i.e. to replace the sum over discrete final states by an integration over a level 
density. Determination of the nuclear level density poses the same problem as in the 
case of compound reactions. Additionally, spectroscopic factors have to be predicted. 
Previously, these were assumed to be constant [32] or energy-dependent [UH HU 02] . 

The spectroscopic factors describe the overlap between the antisymmetrized wave 
functions of target+nucleon and the final state. The number of final state configurations 
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Figure 2. Sums of the contributions to the total direct capture cross sections with 
the averaged (statistical) DC model and the standard DC for 136 Xe(n,7) as function 
of excitation energy [41] 



increases with increasing excitation energy and the overlap of initial and final state 
wavefunctions decreases. Thus, also the spectroscopic factor S decreases. In a simple 
approach, the energy dependence of the spectroscopic factor can be parameterized by a 
Fermi function with 



and the parameters E*,A. Despite of its simplicity, direct capture cross sections can 
be well described by this approach [JT]. Figure [2] shows a comparison between averaged 
direct neutron capture and standard neutron capture calculated in the Lane-Lynn model 
[H] for 136 Xe and using the spectroscopic factors shown in Figure [3j For one-nucleon 
capture, S can also be derived from the occupation probability v 2 for the target state 
into which the nucleon is captured. It can be calculated, e.g., by employing BCS pi] 
or Lipkin-Nogami pairing [38, [37]. For the case of one-nucleon capture on even-even 
nuclei, the spectroscopic factor S for capture can be reduced to S = 1 — v 2 [38l US]- 

Recently, it was found that thermal excitation of target nuclei reduces the sensitivity 
to the parity dependence in the nuclear level density [I]. We can expect a similar effect 
for direct capture, although the number of possible transitions is more limited [ij. A 
detailed discussion will be presented in a forthcoming paper [46J, along with large-scale 
calculations in the averaged direct capture model with excitation-energy dependent level 
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Figure 3. Comparison between experimental and averaged (equation [T]) spectroscopic 
factors as function of excitation energy for 136 Xe+n [41] 



densities and spectroscopic factors. 



4. Summary and conclusion 

The sensitivities of reaction cross sections to predicted nuclear properties were discussed 
here for compound and direct reactions. Obviously, strong single resonances in the 
relevant energy range will also have considerable impact on the resulting astrophysical 
reaction rates. Because the resonance energies have to be predicted to a few hundred 
keV, this poses a continuous challenge to nuclear theory. For a detailed study of the 
actual impact of rate variations on nucleosynthesis in massive stars for nuclides with 
28 < A < 80, see 
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